Analog Computer Applications 
Application Note #3 
Undergraduate Engineering Physics Laser-Themed STEM Project 





Destination Europa: 
Laser Power-Beaming Simulation 


1 Introduction 


This application note is based on a hypothetical space mission to explore just 
beneath the surface of Europa’s ocean at a fixed depth. A submersible probe is 
designed to experience neutral buoyancy (weight = buoyancy). Thus, 
assuming no external horizontal lateral forces acting on the probe, then only 
thrust and drag need to be considered. For this simulation, thrust is assumed 
to bea linear function of laser power. Since a laser is the only available source 
of power for the probe’s port and starboard electric motors (see figure 1), the 
issue of laser beam transmission through seawater must be considered. The 
laser power delivered to the photovoltaic array located at the aft end of the 
probe is assumed to be an exponential function. 


Initially, the probe is given a mechanical ‘nudge’ (not shown in the figure) 
with a constant speed over a short distance just prior to activating the laser. 


This simulation solves and displays the speed of the probe as it travels 
through this alien ocean! 
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Figure 1: Laser Power-Beaming (rendition) 


2 Mathematical modeling 


Starting with Newton’s second law of motion, 


Max = )iforcess-direction = -drag + thrust = -D + T (1) 
May = )iforcesy-direction = 0 (no lateral forces present) (2) 
Ma, = )iforces;-direction = -weight + buoyancy = 0 (neutral buoyancy) (3) 


Returning to (1), 
Max = -D + T 
Assuming 
D= ¥%2CppAv? = «kv? (quadratic drag model), 
T = cP =cPmmaxe*(-OX) = Tmax e*(-ax) = te*(-ax), and 
ax = dv/dt = dx/dt dv/dx = vv’ (per chain rule for differentiation), then 
Mvv’ = -«v? + te’ (-ax) or 
Vv = (-«/Mv + t/Me*’(-ax)/v) with vo = v(Xo) (4) 
where 
M = mass of probe 
x = distance between laser aperture and probe’s photovoltaic array 
v = speed of probe 
« = drag factor 
T= maximum thrust 
a = absorption coefficient 


After several trial function ‘guesses’, the analytical solution takes the form 


v =yv(x)e*(-ax/2) (5) 


Differentiating (5) to obtain the left-hand side of (4), 
v= V(x) (C-a/2 Je®-ax/2)) + e* (-ax/2)(y/(2V (x) 
Equating this expression to the right-hand side of (4), 
yV (x) (-a./2)e* (-ax/2)) + e* (-ax/2)(y/(2V(x)) = 
(-«/M)yvV(x)e* (-ax/2) + (t/M)e*(-ax/2)/@vV(x) 
After a bit of algebra, 
(-./2)x + (1/2) = (-«/M)x + t/(My?) 
Comparing coefficients of like terms, 
-a./2 = -«/M or Ma. = 2« and 
1/2 = t/(My?) or y = V(2t/M) 
Thus, in general, 
v= yv(x)e’(-ax/2) with v(xo) = vo and y = V(2t/M) (6) 


It is useful, and perhaps instructive, to analytically determine the distance at 
which the probe acquires maximum speed and what that speed is. 


Differentiating (6) with respect to x, 
v= V(x)C-a/2)e*(-ax/2) + ye*(-ax/2)/(2V(x)) 
Setting v’ = 0 to determine x*, 

0 = yV(x*)(-a/2)e*(-ax*/2) + ye* (-ax*/2)/(2V(x*)) 
0 =-av(x*) + 1/V(x*) 


x* = 1/oa and 
Vmax = V(2t/M)V(x*)e4 (-ax*) 

= V(2tx*/M)e* (-ax*/2) 

= V(2t/(Ma))e*(-1/2) 

= V(2t/(2k))e*(-1/2) = V(t/«)e*(-1/2) 
Letting a = 0.0641 m1,t = 1.00 N, « = 0.0400 kg/m, and M = 1.247 kg, 
x* = 1/0.0641m1 = 15.6 m, and 
Vinax = V(1.00 N/0.0400 kg/m)e*(-1/2) = 3.03 m/s 
For the computer setup (figure 2), rewrite (4) as 
Vv’ = (-Kv + te*(-ax)/v)(1/M) with vo = v(xo) 

Parameter calculations 
Substituting values into (6), 
vo = v(0.5 m) = V(2 x 1.00 N/1.247 kg x 0.5 m)e*(-0.0641 m™ x 0.5 m/2) 
= 0.881 m/s and 


To = T(0.5 m) = 1.00 N e4(-0.0641 m*™ x 0.5 m) = 0.968 N 


3 Computer setup 
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Figure 2: Computer setup for laser power-beaming simulation 























Parameter | Value 
K 0.0400 
To 0.968 
a 0.0641 
Vo 0.881 
1/M 0.802 











Table 1: Parameter settings for laser power-beaming simulation 


4 Results 





Figure 3: Probe speed as a function of distance traveled* 





Analog Computation Analog Computation Analytical 
x (mm) v (m/s) v (m/s) 
Oscilloscope Oscilloscope Calculator 


Scale: 1m per small division | Scale: 0.2 m/s per small division | Rounded to two 
decimal places 
00.5 0.90 0.88 = vo 

15.6 3.15 3.03 = Vmax 




















Table 2: Solution Comparison: Analog Computation vs Analytical 


*For this application note, the display was produced during a single run by a 
differential equation analog computer prototype using discrete components 
with tolerances between 1% and 10%. 
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